Robustness of spin-triplet paring and the singlet-triplet paring-crossover in 

superconductor / ferromagnet hybrids 
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We have investigated the proximity effect in superconductor/ferromagnet junctions in a systematic 
manner to discuss the relation between the zero energy peak (ZEP) of the local density of states 
(LDOS) and spin-triplet odd-frequency paring. By exactly solving the nonlinear Usadel equations, 
we have found that the ZEP is realized in a wide range of geometrical and material parameters in the 
case of the noncoUinear magnetization. This strongly suggests the robustness of the ZEP induced 
by spin-triplet odd-frequency paring in such systems. We also found that the crossover from singlet 
to triplet paring can be detected by measuring the F layer thickness dependence of the ZEP height. 
Further, we show how to observe signatures of spin-triplet odd-frequency paring and the paring 
crossover by the LDOS measurements. Our results provide a direct way to experimentally detect 
signatures of odd-frequency paring state. 

PACS numbers: 74.45. -|-c, 74.78. Na, 72.25.-b, 85.75.-d 



I. INTRODUCTION 

The study of proximity effect in superconductor (S) 
- ferromagnet (F) hybrid structures has its long history 
following the first theoretical proposal of the so-called 
TT state in a mesoscopic ring containing SFS Joseph- 
son junction. The penetration depth of Cooper pairs 
into a diffusive normal metal (N) is characterized by the 
length scale — \/hD/2TiT, whereas in a ferromagnet 
this length scale is considerably shorter and is given by 

— \/ /2Ecyi. Here T is temperature, D is the diffu- 
sion constant, and i?cx is the magnitude of the exchange 
potential in the ferromagnet. Since the exchange field af- 
fects differently electrons with opposite spins, spin-singlet 
Cooper pairs are fragile under the exchange potential. In 
addition to a short penetration length, the pairing func- 
tion of spin-singlet pair spatially oscillates with changing 
its sign under the exchange potentialj^i^ which enables 
the TT states in SFS junctions.— Although the tt state 
was predicted theoretically in 1970's, it has been con- 
firmed experimentally only recently^"— The details of re- 
search progress have been summarized in several review 
papers i^"— 

Bergeret, Volkov and Efeto\^ proposed theoretically 
new type of proximity- induced superconducting state in 
ferromagnets, the so called long-range spin-triplet pair- 
ing state. Inhomogeneous magnetic structures near the 
SF interface (see Fig. 1) rotate the spin direction of an 
electron, which induces equal spin-triplet s-wave Cooper 



change potential, they have the long range penetration 
length into ferromagnets characterized by ^t- Experi- 
mentally such effect has been observed first as the long- 
range Josephson coupling in SFS junctions;^ where the 
ferromagnet is a half-metallic Cr02 compound. In clean 
half-metallic SFS junctions the equal-spin triplet pairs 
can have the odd-parity even-frequency symmetry^ 
whereas in the experimentally relevant dirty limit case 
the equal-spin triplet even-parity s-wave odd-frequency 
symmetry is the only possible choice. The experimenl^i^ 
as well as anomalous conductance oscillations observed 
in SF hybrids^! stimulated much theoretica l^^d^i^° and 
experimental work^i^— . As a result, a number of signa- 
tures of triplet correlations have been observed. However, 
no experiment so far has reported direct and unequivocal 
evidences of long-range spin-triplet odd-frequency pairs. 

Existing theoriesiiiii"— predict that the presence of 
odd-frequency pairs causes the enhancement of the zero- 
energy local density of states (LDOS)j ^^d''' '— According 



DOS probed 



x = x = Lf 



X — Lp 



pairs in ferromagnets 
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Such Cooper pairs should 



have the odd-frequency symmetry to satisfy the require- 
ment from the Fermi-Dirac statistics of electrons. Since 
equal spin-triplet pairs are not suppressed by the ex- 



FIG. 1: (Color online) Model of an SF junction for observa- 
tion of spin-triplet paring through the local density of states 
(LDOS) measurement. The magnetization directions in Fl 
and F2 layers are coUinear {a = 0) or noncoUinear (a = ■7r/2). 
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to a number of theoretical papers on the proximity effect 
in various SF hybrid structures)^"— the relative fraction 
of odd-frequency pairs to even-frequency pairs depends 
sensitively on junction parameters such as the resistivity 
of F, the transparency of SF interface, the amplitude of 
the exchange energy £^cx, and the geometry of junctions. 
To observe clear evidence of spin-triplet odd-frequency 
pairs in experiments, theoretical studies should show a 
way how to optimize the fraction of odd-frequency spin- 
triplet pairs in wide parameter range tunable in actual 
experiments. 

In this paper, by solving full-spin Usadel equation^i^ 
in a wide parameter range, we systematically calculate 
the LDOS at a surface of a diffusive ferromagnet con- 
nected with a metallic superconductor. In particular, 
we focus on a relation between the magnitude of the 
zero-energy peak (ZEP) in the LDOS and the fraction 
of triplet odd-frequency pairs and show the robustness of 
the presence of the ZEP. It should be noted that ZEP for- 
mation has been reported in measured tunneling conduc- 
tance spectra in oxide-based SF hetero-structures with 
nonuniform ferromagnetsi^"— However the physical ori- 
gin of ZEP in these structures is still unclear. Therefore 
we propose an experimental method to explicitly detect 
a signature of spin-triplet paring by measuring ZEP. 

This paper is organized as follows. In Sec. II, we 
present a model of an SF junction and describe the nu- 
merical method to solve the nonlinear Usadel equation 
in such system. The numerical results of the LDOS for 
various parameters and discussion of the robustness of 
ZEP are presented in Sec. III. In Sec. IV, the summary 
of our results is presented. Throughout the paper we 
confine ourselves to the regime of zero temperature and 
put ks = 1- 



II. NONLINEAR USADEL EQUATION 

In this section, we will introduce a model of an SF 
junction and describe how to numerically solve the non- 
linear Usadel equation with inhomogeneous magnetiza- 
tion in a ferromagnet F. The system we consider is an 
one-dimensional diffusive SF junction where magnetiza- 
tion in the F layer is either homogeneous or rotating like 
illustrated in Fig. 1. The rotation is introduced by di- 
viding the F layer with thickness Lp into Fl and F2 
layers where magnetization are in general non-coUinear. 
The thicknesses of Fl and F2 layers, the SF barrier resis- 
tance, and the F layer resistance are respectively defined 
by Lpi, Lp2, Rb, and i?jv. The exchange field in the F2 
layer is aligned along the z axis, while that in the Fl layer 
is misaligned by angle a in the x-z plane. Importantly, 
the misorientation angle a in exchange-spring ferromag- 



net 



,61-63 



can be experimentally controlled by applying 
an in-plane magnetic field. In our calculation we have 
assumed that the exchange energy i?ex in the Fl and 
F2 layers are the same, and the Fl and F2 interface is 
transparent. 



Because of the spin-rotation of electrons in ferromag- 
nets, we have to solve 4x4 matrix Usadel equation^ 
given by. 



ihD- 



dx 



9-r9 

ax 



\H,g\ -0, 



(1) 



where D is the diffusion constant in F. Here Hamiltonian 
H and the Green function g are respectively defined by 



H = 



g{x,E) 



E{x,E) A{x,E) 
aIx,E) E{x,E) 

g{x,E) f{x,E) 
f{x,E) g{x,E) 



where 



E{x,E) =Ei - Vix) ■ 
E{x,E)^{Eix,-E)y , 

A{x,E) =|a(x,-£;)|* . 



(2) 
(3) 

(4) 
(5) 

(6) 



Here dj with j ~ 1, 2, 3 are Pauli matrices and (Jq = 1 is 
the 2x2 unit matrix. The magnetic moment V{x) in a 
ferromagnet is defined as 

, _ J i?ox(sin a, 0, cos a) for 0<x< Lpi 
^ ~ \ Ec^lo, 0, 1) for Lfi < X < Lp ' ^ ' 

Throughout this paper, • " • and • " • indicate 4x4 and 2x2 
matrices, respectively. In what follows, we only consider 
the s-wave spin-singlet pair potential in a superconduc- 
tor, (i.e., A = Aoi(T2). The particle-hole symmetry re- 
sults in 



gix,E) = -{g{x,-E)}* , 
f{x,E)^~[fix,-E)y . 



(8) 
(9) 



To solve Eq. Ill), we use the Riccati parameteriza- 
tion^"— for the Green function, i. e.. 



g{x,E) 



where 
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7V=(i+77)-\ 

iy=(i + 77)-i. 



(10) 

(11) 

(12) 



Here Riccati parameters 7(x, E) and 7(2:, E) are inde- 
pendent from each other. The normalization condition 
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of the Green function is automatically satisfied under 
the parameterization, i. e., gg = 1. The derivative of 
the inverse matrix dxN can be done as follows 



where 



d^N = - NAN, 



(13) 



(14) 



which is obtained from an identity dx{NN^^) = 0. 

Finally, the Usadel equation ([1]) is reduced to two par- 
tial differential equations for 7 and 7, 



ihD 



dll - {dxi) f {dxi) 



ihD 



dll - {dxi) hd.l) 



+ jE + A - 7A7 = 0, 
(15) 

- i;7 + 7^; - A - 7A7 = 0. 

(16) 



After taking the complex conjugation and E —E 
in the above equations, we find that —j*{x,—E) and 
7(x, E) obeys the same equation. Thus we conclude that 



By numerically solving nonlinear differential equations 
([T5|) and ([T5)) together with the boundary conditions ([2^ 
and (f23| . we will calculate local density of state (LDOS) 



(24) 



and the pair function matrix 



0-2 



fiE)^[fo{E)i + f^ 

ifi{E)+f2{E) 
-ih{E)+iUE) 



-ih{E) - 
-^h{E) 



^fo{E) 
-f2{E) 



(25) 



where Ni^ is the normal DOS and / ~ (/i,/2,/3)- 
Here the components /o, /i, /2, and /a correspond, re- 
spectively, to spin singlet [{WD — lit)) /\/2]i equal-spin 
triplets [|tt) - lU) /\/2], [|n) + lU) /\/2)], and hetcro- 
spin triplet paring [([U) + lit)) 1^2]. 

Since the diffusive limit is considered, the singlet com- 
ponent /o has an even-frequency symmetry while the 
triplet components fi (i = 1, 2, 3) have an odd- frequency 
symmetry. It should be noted that the y component of 
the magnetic moment V(a;) is zero in this paper [see 
Eq. ([7])]. As a consequence the equal-spin triplet com- 
ponent /2 is always zero. 



^{x,E)^^r{x,-E). 



(17) 



At the interface between a ferromagnet and a super- 
conductor, the Kupriyanov-Lukichcv boundary condi- 
tion^ connects Green functions on both sides. 



2TBgdxg = [Gs,g\ , 



where 



Gs = 
9s = 



9 si /sO-2 
fs^2 -9 si 

E 

y/E^ - A2 ' 
iA 

VE^ - A2 ■ 



(18) 

(19) 
(20) 
(21) 



Here Tb = Rb/Rn with R% = Rn{L = ^tJ = 
Rn^tJLf, where ^t, = yJhD/2TTTc [T^ is the supercon- 
ducting transition temperature). In the calculation, we 
have assumed the retarded causality in the Green func- 
tion. By using Eq. (|13|) . we obtain two boundary condi- 
tions 

2rBc),7=2(7,7 + /,(7(T2 7-'72), (22) 
2TBdxl^2gsj + fs{l(J2l~d2). (23) 

Since g,{-E) = g*{E) and f,{~E) = f:[E), 
— 7*(a;, —E) and j{x,E) satisfy the same boundary con- 
dition. 



III. RESULTS 

In this section we study the LDOS and the zero-energy 
peak (ZEP) in LDOS for the cases of the uniform and 
nonuniform magnetization in a systematic manner. In 
order to show the robustness of the ZEP induced by 
spin-triplet odd-frequency paring, we calculate phase dia- 
grams of the ZEP and the pair amplitudes as functions of 
several variables. We also show how to detect signatures 
of the long-range triplet-paring experimentally. 



A. Local density of states and pair functions 

Let us first discuss the LDOS for the uniform mag- 
netization case, i.e., a = 0. In the calculation, we 
set E^^/2-kTc = 0.1 and R%/Rb = 0.2 (a moderate 
proximity-effect regime). We show the F layer thickness 
Lp dependence of the LDOS N{E) at the edge of F2 
layer, namely at x = Lp in Figs. 2(a)^(c). When Lp/^r^ 
is very small or the Thouless energy i^xh is much larger 
than Ecx, i-c, i^xh = hD/L^ ^ Ecx, the mini gap is 
formed due to the proximity effect. In this case we can 
neglect magnetic effects and then the magnitude of the 
mini gap is approximately given by^ 



E' 



Th 



Re 



(26) 



as in the case of diffusive S/N(normal metal) junctions. 
By increasing Lp, i?Th(~' EJ.^) is decreased. Eventu- 
ally, if a resonant condition E^x ~ E^^ is satisfied, the 
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FIG. 2: (Color online) The local density of states (LDOS) N{E) at the edge of the F2 layer i.e., x — Lf, as a function of energy 
E and the F layer thickness Lf for an SF junction with (a) an uniform (a = 0) and (d) a nonuniform (a = ''r/2) magnetization 
with Lfi = O.SCt,, Sex/27rTc = 0.1, and R%/Rb = 0.2. iVo, Ao, and Ct„ are the normal state LDOS, the superconducting 
gap, and the coherence length at T = Tc, respectively. Panels (b) and (e) are magnified figures of (a) and (d) near E = Q, 
respectively. Panels (c) and (f) are magnified figures of (a) and (d) around the small Lf regime, respectively. The arrows 
indicate the resonant condition £cx = Smg, where E-mg, is the mini gap in the case of iJcx = 0. 



mini gap is closed completely [see Figs. 2(a) and (c)]<^ 
As a consequence, the zero-energy peak (ZEP) is devel- 
oped. However, if we increase Lp ov decrease i?Th fur- 
ther, the LDOS profile near E = Q becomes almost flat 
because the mini-gap edges move outwards toward the 
superconducting-gap edge. Therefore, the ZEP can be 
realized only near the resonant condition ii^ox = ^'mg'^'^ 

In order to see more clearly this point, we have plotted 
the zero-energy density of state iV(0) and the amplitude 
of the pair functions |/i(0)| (i = 0, 1,3) as a function of 
Lp in Fig. 3(a). In the uniform-magnetization case, long- 



range triplet components /i and /2 are completely ab- 
sent Moreover near the resonant condition i5ex = ^-mg, 
although the fraction of the short-range triplet compo- 
nent /3(0) is larger than the singlet one /o(0), both com- 
ponents coexist with each other. These observations are 
consistent with previous result J^i^^ 

On the other hand, in the case of a nonuniform mag- 
netization (a = 7r/2), the behavior of iV(0) is drastically 
different from the uniform case due to the appearance of 
long-range spin-triplet paring. The LDOS at the edge of 
F2 layer as a function oi Lp is shown in Figs. 2(d)~(f). 
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FIG. 3: (Color online) (a) The zero-energy LDOS N{0) and the absolute value of the zero-energy pair- functions fi{0) as a 
function of the F layer thickness Lf at the edge of the F2 layer {x = Lf)- fo (black dashed line), /s (blue dotted line), and 
/i (red solid line) are the shot-range singlet, the short-range triplet, and the long-range triplet component, respectively. The 
vertical dotted-line is corresponding to the resonant condition Box ~ Eg. The parameters are Lfi ~ O.S^Tc Ec^/2tiTc — 0.1, 
and R%/Rb = 0.2. 



As clearly scon from Figs. 2(d) and (c), not only near the 
resonant condition but also a wide range of Lf , the ZEP 
is developed. 

In order to understand the physical origin of the ZEP, 
we have plotted the zero-energy LDOS iV(0) and the pair 
functions /^(O) at the edge of F2 layer {x = Lp) as a 
function of Lp in Fig. 3(b). If the long-range triplet 
component /i(0) is largely developed in comparison with 
short-range components (/o and /a), the ZEP is highly 
enhanced. Although short range components show an 
exponentially-damped oscillation as a function of Lf, the 
long-range one decays considerably slowly with increas- 
ing Lp. This is due to the fact that long-range pairs 
(Iff) , \\,\.)) have zero center-of-mass momentum as in the 
case of singlet-pairs in conventional SN junctions. There- 
fore in the case of the long F2 layer i. e., Lp2/(,Ta ^ 1, 
one can have almost pure long-range triplet component 
near the edge of F2 layer. 

It is important to note that the dependences of A^(0) 
and /i(0) on Lf in the regime Lf/Ctc ^ 1 are closely 
related to each other as demonstrated in Fig. 3(b). This 
can be explained as follows. If the long-range compo- 
nent /i(0) is fully developed, the zero-energy LDOS is 
approximately given as 

^ = ^Tr [Rc5(0)] « ^/l^imf, (27) 
as can be obtained by using the normalization condition 
f+ff^l- (28) 



Thus the L p dependence of ZEP is closely related to that 
of the long-range triplet component fi{0)J^ Therefore 
we can conclude that a systematic ZEP measurement by 
changing Lp gives a strong evidence for the signature 
of the long-range spin-triplet correlations. In Sec. IV. C 
we will discuss in more detail the ways to discriminate 
experimentally between the short- and long-range triplet 
components. 



B. Robustness of the zero-energy peak 

Heretofore only limiting cases, e.g., a very weak ex- 
change field)^ a very strong exchange field (like half met- 
als) ^^i^ and a very thin F layer thickness (Lf <C CTc)r^ 
have been explored in previous theoretical studies on 
the ZEP in SF junctions. Moreover in order to vali- 
date the linearization approximation for the Usadel equa- 
tion, only the weak proximity- effect regime (equivalently 
Rn/Rb ^ 1) has been explicitly investigated. Therefore 
the natural question to ask is how robust the presence of 
the ZEP induced by spin-triplet odd-frequency paring is in 
reality? To solve this crucial issue, we will calculate the 
zero-energy LDOS A^(0) and the pair amplitudes fi{0) 
by systematically varying several parameters, i. e., (1) 
the exchange energy E^^ and (2) the barrier resistance 
Rb which is a control parameter for the proximity effect 
without the linearized approximation . Then we would 
like to show the robustness for the presence of ZEP in- 
duced by spin-triplet odd-frequency paring. 
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FIG. 4: (Color online) The phase diagram of the zero-energy LDOS A^(0) at x = Lf as a function of Rn/Rb and -Bcx/^Th for an 
SF junction with (a) an uniform (a = 0) and (b) a nonuniform magnetization {a = 7r/2), where Lpi = 0.5^t^ and Lf = 4.0^Tc • 
The mini gap, the zero-energy-peak (ZEP), and the flat phase are respectively defined by the region with N{0)/No < 0.98, 
N{0)/Nq > 1.02, and 1.02 < N{0)/No < 0.98. Panels (c) and (d) show the phase diagram of a most dominant component of 
pair amplitudes at s = Lf for an SF junction with an uniform and a nonuniform magnetization, respectively. Black, blue, 
and red regions are the singlet /o, the short-range triplet /a, and the long-range triplet /3-dominant phase, respectively. The 
dotted lines are corresponding to the resonant condition iJox ~ Emg. 



In Fig. 4 we show the phase diagram of N{Q)/No and 
a most dominant component of pair amplitudes for an 
SF junction with the uniform [Figs. 4(a) and (b)] and 
the nonuniform magnetization-configuration [Figs. 4(c) 
and (d)] as a function of Rn/Rb and Ecx/Etii- In the 
calculation, we have assumed that Lpi — O.S^Tc and 

In the zero-energy LDOS phase diagram, following 
three phases have been realized, i.e., (i) the mini-gap 
phase with N{0)/No « 0, (ii) the ZEP phase with 
N{0)/No > 1, and (iii) the fiat phase with iV(0)/A^o ~ 1- 
In the calculation we have defined the ZEP phase as 
regions with N{0)/No > 1.02 and the flat phase with 
0.98 < N{0)/No < 1.02 as a matter of practical conve- 
nience. As was already discussed in Sec. III. A, in the 
case of the uniform magnetization, only the short-range 
components /q and /s exist. As a consequence, only in 
the vicinity of the resonant condition the ZEP phase is 
developed [see the dotted line in Figs. 4(a) and (b)]. 

In contrast, if we consider a nonuniform magnetiza- 
tion (a = 7r/2), the ZEP phase is appeared in a wide 



range of the parameter region for Ecx/E^Th and Rn/Rb 
as clearly shown in Fig. 4(c). This can be attributed to 
the appearance of the long-range triplet component /i (0) 
[see Fig. 4(d)]. From above results, one can conclude that 
in the SF junction with nonuniform magnetizations, the 
existence of the ZEP induced by spin-triplet paring is very 
robust irrespective of the device and material parameters 
as long as i?ox/£'mg > 1- Therefore the experimental ob- 
servation of the ZEP gives a unequivocal evidence of the 
spin-triplet paring. This is one of the important findings 
in our work. 



C. Zero-energy peak spectroscopy 

In this section we will investigate the ZEP structure 
in more detail and propose an experimental method to 
explicitly detect a signature of long-range triplet par- 
ing through the ZEP measurement, so to say the ZEP 
spectroscopy. Moreover the possibility for observing the 
singlet to triplet crossover via the ZEP spectroscopy is 
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discussed as well. 

We have calculated the deviation of the zero energy 
LDOS iV(0) from the normal value A^Oi *• e-j I'^'^ol = 
|A^(0)/iVo — 1| at a; = Lf by changing Lp- In the calcu- 
lation we have used Lpi — O.S^Te and Ecx/2'kTc — 0.1 
as in Fig. 3. In Fig. 5 we show a log plot of \5vi^\ as 
a function of in the case of a moderate proximity 
regime {R^/Rg = 0.2) and a strong proximity regime 
{R']^/Rb = 1.0). It should be noted that in the strong 
proximity regime, the well-used linearized Usadel ap- 
proach is not justified at all. 

In the case of an uniform magnetization (a = 0), {Svol 
shows the oscillatory damped behavior. This means that 
Svq changes its sign almost periodically with increasing 
Lp. Such a behavior is consistent with previous the- 
oretical prcdictionsi2ri2i and also with experimental rc- 
sults,2^i26 

In contrast, the behavior of {Si^ol for the case of inho- 
mogeneous magnetization (a = 7r/2) is drastically dif- 
ferent from that for the uniform one. Due to the devel- 
opment of the long-range triplet component /i(0), \Siyo\ 
decay very slowly as a function of i^- And apparently 
the sign of Suq is always positive as long as E'ox > -^'mg 
[see also Fig. 3(b)]. In this case, the typical decay length 

is given by = Re \fhDjib , where b is the imaginary 

part of the energy which accounts for inelastic process, 
for instance, the thermal excitation or the spin-orbit in- 
teraction. Throughout this paper wc set (5/Ao = 10^^. 
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FIG. 5: (Color online) The zero-energy peak spectroscopy. 
The deviation of the zero energy LDOS A'^(O) from the normal 
value A^o, e., Svo = N{0)/No — 1 at as a function of the F 
layer thickness Lp for an SF junction with a = (the black 
dashed line) and a = Tr/2 (the red solid and blue dotted line) 
for different values of R%/Rb- The LDOS is evaluated at 
X = Lp. Parameters are Lpi = 0.5^Tc and _Eex/27rTc = 0.1. 



From above results, we can conclude that by systemat- 
ically measuring the zero energy LDOS as a function of 
Lp,wc can explicitly distinguish between long-range (/i) 
and short-range parings (/o and /a). 

More interestingly, we have found a following peculiar 
feature in the strong proximity regime {R%/Rb = 1.0). 
As clearly seen from the blue dotted curve in Fig. 5, a 
crossover from the short-range (0 < Lp /^p^ ^ 11) to the 
long-range behavior [Lp/S^p^ ^ H) appears. In order to 
see more explicitly this crossover, we have also calculated 
pair functions as a function of Lp in Fig. 6 and found that 
dip positions of ISi^ol Si,tc almost identical with crossover 
points between singlet- and triplet-dominant phases. It 
is important to note that such a singlet-to-triplet paring 
crossover can be regarded as an even-to-odd frequency 
paring one. 

A physical origin of above remarkable phenomena can 
be explained as follows. If all the components fi{i = 
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FIG. 6: (Color online) (a) The deviation of the zero-energy 
LDOS A'^(O) from the normal value A^o, i- e., \Svo\ = 
|A'^(0)/A'^o — 1| at a:: = Lp as a function of the F layer thick- 
ness Lp for an SF junction with a non- uniform magnetiza- 
tion (a — p/2) in a strong proximity regime {R%/Rb = 1.0). 
Panel (b) shows the absolute value of pair functions fi{0) 
at s = Lf as a function of Lp. fo (black dashed line), /a 
(blue dotted line), and /i (red solid line) are the shot-range 
singlet, the short-range triplet, and the long-range triplet 
component, respectively. Parameters are Lpi = 0.5^Tc and 
Ec^/2pTc = 0.1. The gay, blue and red region are correspond- 
ing to singlet-, short-range triplet-, and long-range triplet- 
dominant phases. The upper table shows paring symmetries 
for a most dominant component and the sign of Svq. 
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by assuming \fi \ <^ 1, and using the facts that lm/o(0) = 
for the singlet component and Re/i(0) = for triplet 
components (i = 1, 2, 3) i^°'^^ Therefore, singlet and 
triplet components respectively give a negative and pos- 
itive contribution to Si>o- So if the dominant component 
is changed from singlet to triplet (or from triplet to sin- 
glet) by increasing Lp, Svq changes its sign. This gives 
rise to a dip structure in the l^z^o] vs Lp curve as shown 
in Fig. 6. 

It is important to note that even in the weak or mod- 
erate proximity regime {R%/Rb < 1) such a singlet- 
triplet crossover phenomena can be realized. As was 
demonstrated in Fig. 3(b), the long-range component 
/i(0) is largely enhanced in compared with strong prox- 
imity cases [Fig. 6(b)]. Therefore single crossover (from 
a short-range singlet to short-range triplet) is realized 
only at around the resonant condition. Similarly we 
can regard well known zero-energy LDOS oscillation as 
a function Lp in an uniform SF junction as a multi- 
ple (short-range) singlet to (short-range) triplet crossover 
[see Figs. 3(a) and 5]. 

Therefore we can explicitly conclude that we can iden- 
tify crossover points of the paring symmetry from dip po- 
sitions of I (5 1^0 1 ■ This remarkable feature has never been 
recognized in previous studies and is one of most impor- 
tant findings in our paper. From above results it is clear 
that the systematic LDOS measurement by changing the 
exchange field E^^, the F layer thickness Lp, and the 
barrier resistance Rb gives a clear and unambiguous ev- 
idence of the novel long-range triplet odd-frequency par- 
ing. 

Finally we would like to propose a systematic method 
to detect the signature of odd-frequency paring and 
observe the singlet-to-triplct crossover experimentally. 
Fig. 7(a) shows a scheme of an SF junction in contact 
with an STM tip to measure the differential conduc- 
tance or the LDOS in the F layer. The spatial depen- 
dence of the LDOS of F in the strong proximity regime 
is plotted in Fig. 7(b). As clearly seen from Fig. 7(b), 
the measurement of the position dependence of the zero- 
energy LDOS would enable clear identification of long- 
range odd-frequency paring as well as the the singlet- 
to-triplet crossover. It should be noted that the spatial 
dependence of the LDOS of an inhomogeneous SF junc- 
tions for a weak proximity and small exchange field case 
has been investigated by Cottet4i 



IV. SUMMARY 

To conclude we have systematically investigated the 
superconducting proximity effect in SF junctions with 
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FIG. 7: (Color online) (a) Scheme of an inhomogeneous SF 
junction in contact with an STM tip. Panel (b) shows the 
position X dependence of SvoiO) and pair functions /i(0) 
(j = 0, 1,3) in the F layer, where ± is corresponding to the 
sign of Si/Q. fo (black dashed line), /s (blue dotted line), 
and /i (red solid line) are the shot-range singlet, the short- 
range triplet, and the long-range triplet component, respec- 
tively. Parameters are a — 7r/2, Lfi = O.S^Tc Lf2 = 14.5^Tc, 
R%/Rb = 2.0, and E,^/2-kT^ = 0.1. 



an uniform and a nonuniform magnetization in terms of 
spin-triplet odd-frequency paring. By solving the nonlin- 
ear Usadel equation fully numerically, we have calculated 
the LDOS in a fcrromagnet and found following remark- 
able results. 

(l)In contrast to the case of the uniform magnetiza- 



tion 



70.71 



the LDOS in SF junctions with a nonuniform 
magnetization has a zero energy peak in a wide range 
of parameters, indicating the robust presence of the ZEP 
induced by the spin-triplet odd-frequency paring. 

(2) The ZEP height is damped very slowly with increas- 
ing L p due to the development of long-range spin-triplet 
paring. This behavior is very contrast to uniform magne- 
tization cases in which the zero energy LDOS shows ex- 
ponentially damped-oscillation as a function of Lpi^— 

(3) The dip position of \5vq\ is corresponding to the 
crossover-point between singlet and triplet or even and 
odd paring. This means that the ZEP spectroscopy gives 
us clear information on the paring-symmetry of Cooper 
pairs. 

Above remarkable results clearly indicate that the ex- 
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pcrimcntal observation of the ZEP for SF junctions with 
a non-uniform magnetization provides a smoking gun for 
the existence of the novel spin-triplet odd-frequency par- 
ing. 

In this paper, we have discussed proximity effect as- 
suming spin-singlet s-wave superconductor as a bulk 
state of S. An extension to unconventional superconduc- 
tors is possible based on more general boundary condi- 
tion JliI^ taking account of the Andreev bound state. 
Proximity effect in spin-triplet p-wave superconductors 
is interesting22r^ since odd-frequency pairing is also in- 
duced from bulk superconductor without the exchange 
energyj^ In addition, we have especially focused on the 
LDOS. It is interesting to discuss the anomalous Meiss- 
ner effec t^^i^^ and the surface impedance^ due to the 
proximity effect by odd-frequency pairing. 
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